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In a Bernal-stacked graphene bilayer, an electronic state in Landau level N = is described by its 
guiding-center index X (in the Landau gauge) and by its valley, spin, and orbital indices £ = ±K, a = 
±1, and n = 0, 1. When Coulomb interaction is taken into account, the chiral two-dimensional 
electron gas (C2DEG) in this system can support a variety of quantum Hall ferromagnetic (QHF) 
ground states where the spins and/or valley pseudospins and/or orbital pseudospins collectively 
align in space. In this work, we give a comprehensive account of the phase diagram of the C2DEG 
at integer filling factors v € [—3, 3] in Landau level N — when an electrical potential difference 
Ab between the two layers is varied. We consider states with or without layer, spin, or orbital 
coherence. For each phase, we discuss the behavior of the transport gap as a function of Ab, the 
spectrum of collective excitations and the optical absorption due to orbital pseudospin-wave modes. 
We also study the effect of an external in-plane electric field on a coherent state that has both valley 
and spin coherence and show that it is possible, in such a state, to control the spin polarization by 
varying the strength of the external in-plane electric field. 



PACS numbers: 73.21.-b,73.22.Gk,72.80.Vp 



I. INTRODUCTION 

Electrons in a Bernal-stacked graphene bilayeri^ be- 
have as a chiral two-dimensional gas of massive Dirac 
fermions (C2DEG) 3 . The chiral nature of the electrons 
lead to transport and optical properties that are different 
from those of conventional semiconductor 2DEG'S or of 
the 2DEG in monolayer graphene. In particular, in the 
absence of Coulomb interaction, the Landau level (LL) 
spectrum is given by En — ±Hu}*y/N (N + 1), where 
the Landau level index N = 0, ±1,±2, ... and the effec- 
tive cyclotron frequency lo* — eB/m*c where B is the 
magnetic field and m* the effective mass of the electrons. 
Each Landau level is four-time degenerate when counting 
valley and spin degrees of freedom with the exception of 
Landau level N = which is eight-time degenerate. In- 
deed, an electronic state in Landau level N = is speci- 
fied by its guiding-center X (in the Landau gauge), spin 
a = ±1, valley £ = ±K and orbital n = 0, 1 indices^. 
(In Landau level N = 0, valley and layer degrees of free- 
dom are equivalent.) When the Coulomb interaction is 
negligible with respect to the disorder broadening at low 
temperature and when the small Zeeman splitting is ne- 
glected, the eight states in N = are degenerate and the 
Hall conductivity has plateaus at a xy = 4Me 2 /h where 
M = ±1,±2,..A 

In recent transport experiments 5-12 , it was shown that 
in sufficiently pure sample, when disorder is low or when 
the magnetic field is large enough, the Coulomb inter- 
action completely lifts the degeneracy of the N = 
octet and lead to the formation of seven new plateaus 
in the Hall conductivity i.e. a xy — vMe 2 /h, where 
v G [— 3,3] . These plateaus were attributed to the for- 
mation of broken-symmetry many-body ground states. 
These states can alternatively be described as quantum 
Hall ferromagnets (QHF's) where the spin and/or val- 
ley pseudospins and/or orbital pseudospins are sponta- 



neously and collectively aligned in spaced. 

In bilayer graphene, a top-bottom gates voltage imbal- 
ance can be applied to create a potential difference Ab 
(which we call the "bias" hereafter) between the two lay- 
ers. In dual-gated bilayer graphene, Ab and the total 
density of electrons in the bilayer can be controlled inde- 
pendently. This allows the phase diagram of the C2DEG 
to be studied as a function of Ab, magnetic field and 
temperature. Such study has been done by Weitz et 
alJ^ in high-quality bilayer graphene suspended between 
a top gate electrode and the substrate. The measure- 
ments show a series of phase transitions between different 
QHF states as As is increased at given filling factor and 
magnetic field. Special attention has been given to the 
filling factor v — where the precise nature of the ground 
state when B —> near zero bias is still debated-^. All 
experiments were done at relatively small magnetic field 
B < 10 T with the exception of the experiments reported 
m Ref. | where B reached 35 T. 

Various aspects of the QHF states in bilayer graphene 
(in particular the nature and the evolution of the 
ground state of the C2DEG near charge neutrality as 
B — > 0) have been studied theoretically by a number 
of authors^—. In the work of Gorbar et al.— and 
Shizuya^, the phase diagram of the C2DEG as a func- 
tion of bias for v = 0, 1, 2, 3 is presented. Gorbar et al. 
have considered the effect of both statioii and dynamical 
screenings^ of the Coulomb interaction. The modified 
gap equation captures the linear scaling of the transport 
gaps with the magnetic field which is seen in all the exper- 
iments at low magnetic field. The transport gaps are also 
strongly reduced with respect to the unscreened case and 
become more comparable to those observed experimen- 
tally. The dynamical screening was found to reproduce 
the offset in the behavior of the gap with magnetic field 
in the spin-polarized QHF state near zero bias which is 
seen in the experiments^. Shizuya points out that the 
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filled levels N < — 1 cannot be considered as completely 
inert. Instead, they lead to a correction A„ of the energy 
of the orbital levels n — 0, 1 that can change the ordering 
of these states in a way that depend on their occupation. 
A negative capacitance effect is also found that suppress 
rotation of the valley pseudospins. 

The work that we present in this paper extends our 
previous study of the phase diagram of the C2DEG at 
zero bia o 13 i 24 ~ — and complements the study of Gorbar 
et al. and Shizuya. We give a comprehensive account 
of the phase diagram of the C2DEG in bilayer graphene 
at all integer filling factors v G [—3,3] in Landau level 
N = as a function of an applied bias at a fixed 
magnetic field. Our analysis is based on an effective 
two-band model 2 ^ which describes the low-energy physics 
near the valleys K±. We explicitly take into account as 
symmetry-breaking terms the Zeeman splitting Az and 
the bias Ab- In the effective two-band model, the bias 
As lifts the degeneracy between the orbital levels n = 
and n = 1 by a small amount £,/3Ab where j3 = 
with 71 the interlayer hopping between carbon atoms 
that are part of a dimer. We also include in our model 
the interlayer next nearest-neighbor hopping term 74 be- 
tween carbons atoms in the same sublattices. This term 
causes a small asymmetry in the electronic band struc- 
ture and was neglected in previous studie o 17 i 19 . It com- 
bines with the correction £/3Ab to give an energy differ- 
ence w 2/37174/70 + £(3Ab between the n = 1 and n = 
orbital levels where 70 is the intralayer hopping energy 
between nearest-neighbors. This correction is thus finite 
at zero bias and breaks the orbital degeneracy. Our phase 
diagram is not electron- hole symmetric around v = and 
the sequence of phase transitions is different for each fill- 
ing factor. 

In our analysis, we treat the electron interaction in 
the Hartree-Fock approximation (HFA) and compute the 
collective excitations and electromagnetic absorption of 
the different phases of the C2DEG in the generalized 
random-phase approximation (GRPA). We include in our 
study both uniform and non-uniform states and allow for 
the possibility of any type of coherent (or QHF) state. By 
coherent state, we mean a state where the average value 

c \,o,n,X c £' >' >' \x) 7^ for £ £' and/or a ^ a' and/or 

n 7^ n! where ct a n x creates an electron in a state with 
quantum numbers £, cr, n, X. In our phase diagram which 
is summarized on Fig. 5, the layer- coherent states occur 
at very small bias because of the small interlayer dis- 
tance d = 0.34 nm in bilayer graphene. As the bias is 
increased, we find around a critical bias corresponding 
to the regions where the Hall conductivity ceases to be 
quantized in the experiments^ a state with both layer 
and spin coherence. The orbital coherent states occur at 
a much larger bias corresponding to the situation where 
level n = 1 gets lower in energy than level n — in valley 
(see Fig. 0|. In-between these coherent states are 
various incoherent states, some of which have been stud- 
ied before^. Interestingly, we find that the application 



of an electric field in the plane of the layers can produce 
a new state where all three coherences (layer, orbital and 
spin) are present. In such a state, it is possible to control 
the degree of spin polarization by changing the strength 
of the external in-plane electric field. 

We also present a study the properties of the different 
ground states in the phase diagram. For all filling fac- 
tors, we show how the transport gaps evolve with bias. In 
most cases, this evolution is qualitatively similar to that 
obtained with screening corrections^!. We compute all 
the intra-LL collective excitations in the various phases 
showing that all coherent states but the orbital state 
are characterized by a linearly-dispersing gapless (in the 
long-wavelength limit) Goldstone mode. This mode be- 
comes gapped after the transition to an adjacent inco- 
herent state. In the orbital phase, the orbital-pseudospin 
Goldstone mode dispersion is anisotropic and this mode 
becomes unstable at a finite wave vector indicating a 
transition to a charge-density-wave stat o 25 ' 26 . We iden- 
tify the number of spin-waves and orbital modes in each 
phase. These later modes are active in optical absorp- 
tion. The inter-LL and some intra-LL magnetoexcitons 
have been computed recently 2 ^— and we comment on 
the difference with our results and how the presence of 
the inter-LL magnetoexcitons in the spectrum may com- 
plicate the detection of the intra-LL excitations. The 
main results of our paper are summarized in Fig. 5 (phase 
diagram), Fig. 8 (transport gaps) and Fig. 9,10 (collec- 
tive mode dispersions). 

This paper is organized in the following way. Section 
II introduces the two-band model of bilayer graphene 
with the resulting LL spectrum in finite magnetic field. 
Section III summarizes the Hartree-Fock and generalized 
random-phase approximations that we use to take into 
account the Coulomb interaction and gives the formal- 
ism for the calculation of the electromagnetic absorption. 
Our numerical results for the phase diagram, transport 
gaps, collective excitations and optical absorption are 
presented in Sec. IV. In Sec. V, we show how the ap- 
plication of an in-plane electric field allow to control the 
spin polarization in some phases. We conclude in Sec. 
VI. 



II. TWO-BAND MODEL OF BILAYER 
GRAPHENE 

A. Crystal structure and tight-binding 
Hamiltonian 

The crystal structure of a Bernal-stacked graphene bi- 
layer is shown in Fig. [TJ Each graphene layer is a two- 
dimensional crystal with a honeycomb lattice structure. 
The honeycomb lattice can be described as a triangu- 
lar Bravais lattice with a basis of two carbon atoms 
A n and B n where n = 1,2 is the layer index. The 
two basis vectors are given by ai = ao (1/2, — v3/2)and 
a 2 = cto (1)0) , where ao = 2.46 A= y/3c is the lattice 
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constant of the underlying triangular Bravais lattice and 
c = 1.42 A is the separation between two adjacent carbon 
atoms. The distance between the two graphene layers is 
d = 3.4 A. In the Bernal stacking arrangement, the upper 
A sublattice is directly on top of the lower B sublattice 
while the upper B sublattice is above the center of a 
hexagonal plaquette of the lower layer. 




hopping between carbon atoms that are part of a dimcr 
(i.e. A\ — B 2 , these sites are called the high-energy sites), 
73 the interlayer NN hopping term between carbon atoms 
of different sublattices (i.e. A 2 — B\) and 74 the inter- 
layer next NN hopping term between carbons atoms in 
the same sublattice (i.e. A\ — A 2 and B\ — B 2 ). The 
energy S represents the difference in the crystal field be- 
tween sites Ai,B 2 and sites A 2 ,B\. In this work, we 
neglect the trigonal warping term 73, a correct approxi- 
mation at sufficiently high magnetic fields. 
If we define the spinor 
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l,k,' 



l 2.W., 



(2) 



where a] 



creates an electron on site A{B) in 



layer i = 1, 2 with wave vector k and spin a = ±1, then 
we can write the second-quantized tight-binding Hamil- 
tonian in the basis (Ai, Bi, A 2 , B 2 ) as 



(3) 



FIG. 1: (Color online) Crystal structure of a Bernal-stacked 
graphene bilayer. 

The Brillouin zone of the reciprocal Bravais lattice is 
shown in Fig. [5] We choose the two nonequivalent valley 
points to be 



Kt = 



where £ = ± is the valley index. 
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FIG. 2: Brillouin zone of the triangular Bravais lattice and 
definition of the valleys indices K±. 

The electronic dispersion is obtained from a tight- 
binding model with the parameters^: 70 the nearest- 
neighbor (NN) hopping in each layer, 71 the interlayer 



with the matrix 

H°Ak) = 

( ±A B + <5 
-\oA z 

-7oA* (k) 



(4) 



\A B 



V 



-74A (k) 
-71 



- 7o A (k) 


-74A (k) 



-74A* (k) 


-\A B 
-\oA z 

~7oA* (k) 



-7i 
-74A* (k) 

-7oA(k) 
\A B + 5 

V 



where Az = g^ B B with g = 2 is the Zeeman energy. 
We have also included in H ° an external transverse elec- 
tric field that creates an electrical potential difference (or 
bias) A B between the layers. The function A (k) is de- 
fined by 



A(k) 



i=i 



(5) 



where the summation is over the vectors connecting a 
site A\ to its three nearest-neighbors in the same plane 
i.e. <5i = a (1/2, 1/2^3) , 8 2 = a (-1/2, 1/2 >/3) , S 3 = 

a (0,-1/ V3)- 

Using A (K^ + k) m — \/3ao (£,k x + ik y ) /2 and setting 
74 = 6 = Ab = Az = 0, we find for the the band 
structure near the points K$ the bands 



Ei,± (p) 
E 2 ,± (p) 



±71 ± 



2m* 



± 



P 



"2m* 



(6) 
(7) 



4 



where the momentum p is measured with respect to fiK^ 
and the effective mass is defined by 



(8) 



The band structure consists of four bands. In the ab- 
sence of bias, the two middle bands meet at the six val- 
ley points. The two high-energy bands are separated by 
a gap 71 from the two middle, low-energy bands. The 
bands Ei,± (p) remain degenerate at p = when 73,74 
and A are finite if Ab = 0. This degeneracy is lifted by 
a finite A^. 

For a neutral bilayer, the chemical potential is at the 
energy E = 0. The low-energy excitations (E « 71) of 
the tight-binding model can be studied using an effective 
two-band mode l 27 ' 28 . This model gives for each valley 



(9) 



V 



1 

2m" 



+ V-HP-P+ 
-\oA z 

{P X + iPy) 



1 

2m ' 



(Px -iPyY 



VtP+P- 
z 



where we used the basis (A2,Bi) for K- and (-81,^2) 
for K + and defined p± = p x ± ip y and 



1 f A B ^ ls S 
2m* V 71 70 71 



(10) 



—i^/nip n -i (x) where ip n {x) with n — 0,1,2,... are the 
eigenfunctions of the one-dimensional harmonic oscilla- 
tor. 

For all calculations done in this paper, we choose 2 - for 
the value of the parameters 



7o 
7i 

74 
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3.1 eV, 
0.39 eV, 
0.12 eV, 
0.0156 eV. 



(16) 
(17) 
(18) 
(19) 



We have checked that the band dispersion obtained with 
this choice of signs for the hopping terms is consistent 
with that reported in the literature™ With the magnetic 
field in Tesla, we have 



P = 8.86 x 10~ 3 B, 
£ = 0.4 04B meV, 



while 



fku* c 
A z 



3. 465 meV, 
0.11585 meV, 



e 



11.25VB meV. 



(20) 
(21) 



(22) 
(23) 

(24) 



In the calculation of a, we take n ss 5 for the effective 
dielectric constant at the position of the graphene layers. 
At B = 10 T, 



In this model, the presence of a quantizing perpen- 
dicular magnetic field is accounted for by making the 
Peierls substitution p — > P = p + eA/c (with e > 0), 
where V x A = B =Bz. Defining the ladder operators 
a = (P x - iPy)£/V2fi and a) = (P x + iP y )£/s/2H with 
the magnetic length £ = ^hc/eB, we get 



where 



and 



tv 



C 2 a 2 
-\oA z 



Ci 
Ci,± 

C2 



P 2™±+8 
V 7o 

Ci±^A B , 



/?7i 1 + 2 




7i 



(11) 



(12) 
(13) 
(14) 

(15) 



The effective cyclotron frequency is lo* = e,B/m*c. 
In Eq. (fTTj) . the ladder operators are defined such 
that a'(p n (x) — i\Jn J r lip n +i (%) and aip n (x) = 



d/a = 3.5911 x lO- 2 y/B = 0.114, (25) 
A z /a = 1.029 1 x 10- 2 VB = 0.0325, (26) 
Hu*/a = 0.307 56VB = 0.973. (27) 



B. Landau levels and eigenstates of the 
non-interacting Hamiltonian 



When 74 = S = As = Az 
spectrum of H® a is given by 



0, the Landau level 



E% = sgn (N) y/\N\ (\N\ + \)hw* 



(28) 



where N = 0, ±1, ±2, ... is the Landau level index and sgn 
is the signum function. The corresponding eigenvectors 
of a given spin are 



1 ( fr|JV|-l,X ( r ) 

y/2 V -sgn( N )h\N\+i,x (r) 



(29) 



for N ^ 0. We use the Landau gauge A = (0, Bx,0) 
where the eigenstates are 



(r) 



1 



(30) 



with X the guiding-center index. All Landau levels 
N =/= are four-fold degenerate including spin and valley 
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degrees of freedom in addition to the guiding-center de- 
generacy AL = S/2ir£ 2 where S is the area of the 2DEG. 
The Landau level N — is an exception because there 
are two degenerate spinors with zero energy which are 
given, in the basis (A 2 ,B 1 ) for K_ and {B\, A 2 ) for K + , 
by 



W 
ho,x (r))'l h hX (r) 



(31) 



It follows that N = is eight-fold degenerate. In this 
paper, we restrict the Hilbert space to the Landau level 
N = and use the index n = 0, 1 to refer to the two 
"orbitals" ip n =o (x) and (f n =i (x). With finite values of 
74, 5, Az or Ab, the valley, spin, and orbital degeneracies 
are lifted and the noninteracting energies become 



- C/ f,<7,n=0 



1 



1 



2 £A B - -aA z , (32) 
^Afl-iaAjf + eSAfl + Ci- (33) 



The corresponding eigenspinors are still given by Eq. 
(|31[). Note that the structure of the sublattice spinors 
in Eqs. (f3"Tj) is such that states from different valleys 
(which are localized on different layers) have no overlap. 
For N = 0, the layer index is thus equivalent to the valley 
index. 



C. Limit of validity of the two-band model 

Fig. [3] shows a comparison between the four-band and 
two-band models for the electronic dispersion in Landau 
levels N = —2,-1,0,1,2 and valley K_ using the val- 
ues of the parameters given previously. The agreement 
between the two models is excellent for N = where 
the difference in energy is of the order of 1%. For levels 
\N\ > 0, the difference in energy between the two models 
is much more important. Note that the two sub-Landau 
levels of N = intersect level N = 1 at A B 0.15 eV. 
For the valley K + (not shown in the figure) the crossing 
occurs at a smaller bias Ab ~ 0.10 eV corresponding to 
an electric field E± sa 300 meV/nm between the layers. 
In our calculation we must keep the bias smaller than 
« 0.10 eV (i.e. As/a < 2.8 for k = 5) for our model to 
be valid. 

Fig. @] shows the ordering of the four levels of a given 
spin in N — at finite bias. The correction £f opens a 
gap between the two orbital states n = and n = 1 which 
is independent of the bias. The effective two-band model 
introduces a correction (3 A b to this gap that has different 
signs in the two valleys as indicated in the figure. When 
combined with £i, the gap in valley K + is positive at 
all biases while the gap in valley K _ changes sign (level 
n = 1 gets below level n = 0) when (3Ab > Ci i- e - f° r 
A B > 0.046 eV (i.e. A B /a = 1.3 for k = 5). 
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FIG. 3: (Color online) Comparison between the complete 
four-band model (symbols) and the approximate two-band 
model (lines) for the electronic dispersion in Landau levels 
N = -2,-1,0,1,2 and valley K-. 



K_,n=l 



K + ,n=l 



K_,n=0 



K + ,n=0 



FIG. 4: Ordering of the four levels of a given spin in Landau 
level N = at finite bias As. 



III. INTERACTING CHIRAL 
TWO-DIMENSIONAL ELECTRON GAS 

We now add the Coulomb interaction to the noninter- 
acting Hamiltonian. Hereafter, we use the same basis 
(A2,Bi) for both valleys and define the field operators 
(r) by 



(r) 



^ ( K, x (r) 



(34) 
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and 



(r) = J2 



x 



h n ,x (r) 




(g) \a) c +!<T>n ,x- (35) 



The second-quantized noninteracting part of the Hamil- 
tonian is given by 

Ho = J2 /**L»( r )^w( r ) ( 36 ) 

= E E ^°,o-,n C e,o-,n,X c e>o-,n,^- 

For the second-quantized Coulomb interaction, 

y = \ E EE/^/^U.w ( 37 ) 

ni,...,n 4 ff,o-' ^ ^ 

where the Coulomb potential 

„2 



^,e' ( r ) 



« |r - r' + (1 - dz| 
has the Fourier transform 
1 x-^ 27re 2 



(38) 



V (r ) ^r^e^'Me^ 1 ^^), (39) 
6 ' Kg 



where q is a two-dimensional vector in the plane of the 
bilayer. The terms that do not conserve the valley index 
in Eq. (|3"T|) are very small and usually neglected^ .. 



A. Hartree-Fock Hamiltonian 

In order to describe the different phases of the uniform 
C2DEG, we define the operators 



a,a' _ 1 \ ~> f 

Pn,n' ~ aj 2-^i C u,n,X Ca ' 
"<p x 



.X-: 



(40) 



where n x {c a , n ,x ) creates (destroys) an electron in a 
state (a,n,X). The index a combines the spin and val- 
ley indices and we use and a a to refer to the specific 
spin or valley index of a = (a^,a a ). The set of aver- 
age values {y°n'n'/} gi yes a complete description of an 
uniform ground state. They are the order parameters 
of that state. The diagonal elements {(pn%)} are the 
filling factors of levels (a, n) , while the off-diagonal ele- 
ments are the "coherences". For nonuniform states, it is 

necessary to define the order parameters |^p°'°/ (G)^| 
where G is a reciprocal lattice vector and (G)^ 

the Fourier transform of (p^^ (r)Y We refer the reader 

to Refs. I25II26I where the formalism for this case is dis- 
cussed in more details. The Hartree-Fock Hamiltonian 



can be written in terms of these operators by (we adopt 
the convention that repeated indices are summed over) 



H HF = N v E a . nP a n % (41) 

-AT Tr( a «' fc s) fnWA* \ r, b ' a 

-l\ v Si. ni ,n 4 ,n 3 ,n 2 [V ) \P ni ,n 2 I Pn 3 ,n 4 ' 



where 



E,,, = K„ +n[u at -^\j, (42) 



with E® n given by Eqs. (|32ll33j) . The Fock interaction is 
defined by 



^ni ,n 2 ,ri3 ,n 4 



(q) 



dpi 2 1 „ 

2tt yl nun2 



(P) 



(43) 



(-P)e 



iqxp^ 2 -pd(l-5 a b ) 



with the form factors 
Kb,o(q) = 
(q) = 

#1,0 (q) = 

Ko.i (q) = 



q 2 f 



' V2 



(44) 
(45) 

(46) 

(47) 



These form factors capture the character of the two dif- 
ferent orbital states. In Eq. (j42|) . v = v + 4 is the 
number of filled levels in N — 0. We reserve the sym- 
bol v e [-3, 3] for the filling factor of the C2DEG. In 
deriving Eq. (|41[) , we have taken into account a neutral- 
izing positive background so that the only contribution 
from the Hartree and background terms is the capacitive 
energy given by the term in parenthesis in Eq. (|42l) . In 
this term, u ai = J2 n .a (Pn i ,n'°' i '' T (0)) is the total filling 
factor in valley aj. Detailed expressions for the Hartree 
(see next section) and Fock interactions H and X are 
given in Appendix A of Ref. I25I . 

The Hartree-Fock energy per electron is given by 

Ehf 



1 



1 d 



~ a,n \rn,n/ 



(48) 



2D 



-X 



n 1 ,n 4 ,n 3 ,n 2 



M 



where N e is the number of electrons in the 2DEG and 
vk ± are the filling factors of the two valleys. 

At q = 0, the only nonzero matrix elements of the Fock 
interactions are 



-^0,0,0,0 (°) 
4,0,0,1 (0) 



Ac,^,f 1 , 1 (0) = i A c , 



*f,to,o (0) 
4,ti,o (0) 



Ac, 
Ac, 



(49) 
(50) 
(51) 
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and the corresponding interlayer terms which must be 
computed numerically. We have defined 



A c = 



(52) 



with 



2 /k£. 



B. Calculation of the order parameters 

We define the single-particle Matsubara Green's func- 
tion 

G a nt n2 (X, t) = - (T T c a ^ x (r) cl n2 X (0)) , (53) 

where T T is the imaginary time ordering operator, such 
that the order parameters are given by 



iV v x 



(54) 



The equation of motion for the Green's function in the 
Hartree-Fock approximation is 



(55) 



where /i is the chemical potential, uj n a fermionic Mat- 
subara frequency and 



jja,c 

"1,13 



,713 ,ni ,ri2 (0)< ) (56) 



are the self-consistent Fock potentials. 

The self-consistent Eq. (1551) can be put in a 8 x 8 ma- 
trix form by defining superindices and then solved nu- 
merically in an iterative way in order to get the order 
parameters. 

The Hartree-Fock equation of motion for the Green's 
function leads to the sum rule (at T = K) 



lp a ' h )\ = (p a ' a ) = v a 

\rn,m/\ \rn,n/ n' 



(57) 



where is the filling factor of the (a, n) level. By defi- 
nition 



(p a ' b ) = (p b ' a y 

\rn,m/ \t m.nf 



(58) 



C. Collective modes in the generalized 
random-phase approximation 

To study the collective excitations, we compute the 
two-particle Green's functions 



Arti ,n2 ,^3 ,ri4 \Mj ' 7 



(59) 



= -N v (T T p a n b un2 (q,r) p^„ 4 (-q, 0)) 

+ ^«,n a (q)><Pn' 3 d >n4 (-q)> 



in the generalized random-phase approximation (GRPA). 
In this approximation, Xn'i'nfn 3 m (q> T ) ls the solution 
of the equation 



^f;„ c ^3,„ 4 (q,^n) (60) 

y (0)o,6,c,d / -q 1 
, £ (0)o,6,e,e / Q ) 

x Hn^ngjnj ,n s (q) Xnf^n^,n 3 ,m (l: 
_I y (0)o,6,e,/ / jq ) 

where f2„ is a bosonic Matsubura frequency and the 
Hartree interaction 



(a ? ,b £ ) 



ni,n2,n3,»4 



(q) = — c Km,n 2 (q)if„ 3 ,„ 4 (-q) (61) 



Xf: 



The two-particle Green's functions xSnt'niiw (Qj^n) 
satisfy the set of equations 



[im n - (E b , n2 - Ea, ni )] (q» «n) (62) 

<Tja,e (0)e,b,c,d ( q ) 



-u: 



e,b (0)a,e : c,d 
?i2,m Xni ,m,n3,n4 



(q, n n ) 



Eq. (|60t can be represented by a set of bubbles (Hartree 
terms) and ladder (Fock terms) diagrams. The function 

Xn°i,n2,n 3 in4 is the Hartree-Fock approximations for the 
two-particle Green's functions. It includes the Hartree- 
Fock self-energy corrections but not the vertex correc- 
tions. Note that two-particle Green's functions depend 
only on the order parameters (Pn b n ) computed in the 
HFA. Eqs. (|6QI62j) can be solved numerically by defining 
superindices and then writing them in a 64 x 64 matrix 
form. The collective excitations are then given by the 
poles of the retarded Green's functions X^,n2,ni%,4, (q> u ) 
which arc obtained by the analytic continuation if2 n — > 
lu + id of the corresponding two-particle Green's func- 
tions. To derive the dispersion relations, we follow these 
poles as the wave vector q is varied. 



D. Pseudospin description 

We showed above that the coherent states of the 
C2DEG can be described by the set of order parame- 
ters l\i°n'n'}} ■ These states are also quantum Hall fer- 
romagnets (QHF's) and can also be described by using a 
pseudospin language where the two valley states £ = + 
(£ = — ) are associated with valley-pseudospin up (down) 
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and the two orbital states n = (n — 1) with orbital- 
pseudospin up (down). 

In this language, the total spin, valley pseudospin, and 
orbital pseudospin components of the electron gas are 
given by 



Si = ^ ( C L,n,X°'a ) ) /9 c €,/3,n,x) , (63) 

L ' = ^rE Y,( c i^x< T [% c i\o l ,n,x) , (64) 



a,n,X 



0i = X]( C L, n ,X CT ni' c «.« J n',x), (65) 

where a^'s are the Pauli matrices and the total filling 
factor is 



(66) 



i,a,n,X 



Note that these 10 fields do not provide a complete de- 
scription of a state. One must also consider the 54 other 
combinations of indices (the 64 order parameters are not 
all independent, however) .We will use both Si, Li, Oi and 
the order parameters (pn b m ) to characterize the ground 
states of the C2DEG. 



E. Induced dipoles 

The coupling of the C2DEG with a uniform external 
electric field in the plane of the layers is given by 



n(r)= J2 CnW^W' ( 68 ) 

cr,£,n,m 

Fourier transforming Eq. (|67p and using the form factors 
defined in Eqs. (|44H47I) . we can show that in an homoge- 
neous stat o 20 i 25 , He gives the dipolar coupling 

H B =-d--E ext , (69) 

with the total dipole vector defined by d = ^ a d a where 

d a = -V2£eN v ( P r (0) x - p a y a (0) y) (70) 

is the dipole moment in valley aj with spin a CT . We have 
defined here 



a 1 / a. a . a,a\ 

Px = 2 0°o,i + Pi,o) 



a, a a,a\ 



2i 



(71) 
(72) 



It is possible to control the orientation of the orbital pseu- 
dospins in the x — y plane with an external electric field. 



F. Electromagnetic absorption 

The total current operator in second quantization is 
given by J = J2 a J awith 



(73) 



The current operator j a (r) is derived from the Hamil- 
tonian in Eq. ^ by making the Peierls substitution 
p — > P = p + eA ext /c and then taking the derivative 
with respect to the external vector potential A. ext 



dH° 



3 a,; 



dAf 



where i = x,y. This gives 



with 



E% = Ci + a^An 



i.e. the bare gap in valley £ and 



(74) 



(75) 



(76) 



(77) 



where H^p = N V E® n p^'n is the noninteracting Hamil- 
tonian. 

To compute the electromagnetic absorption per unit 
area, we define the two-particle current-current Green's 
function 



H E = -e J drn(r)0(r) (67) 
where E e:rf = — V0 (r) . The total density is given by which gives 



XJ a j,(r) = -^(T T J a (r)J p (0)), (78) 



W = (fj^H^b (79) 



x X ^ 6 (q = 0,iO„) 



where a, /3 = x,y and x — y,y — x. Using Eqs. (|71|72j) . 
we have for example 

a,b.c.d ^ / a,b,c,d , a,b,c,d . a,fc,c.d . a,b,c,d 

Xp x ,p x ~ J ^Xl, 0,1,0 Xo, 1,1,0 + Xl, 0,0,1 "+" Xo, 1,0,1 



and similarly for the other components. The absorption 
can only involve these four combinations of orbital in- 
dices whatever the polarization of the electric field of the 
electromagnetic wave. The retarded current-current re- 
sponse function XJ a ,Jf> ( w ) i s obtained from the analytic 
continuation iQ n — > uj + id and the electromagnetic ab- 
sorption for an electromagnetic wave of amplitude Eq 
linearly polarized in the direction a is given by 



P a (u) 



XJ a ,J a H 



(81) 
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This formula is valid at finite frequency only since we 
have neglected the diamagnetic contribution to the cur- 
rent. 



G. Absorption in the incoherent phases 

If there is no coherence in a phase, then 

(pn',m) = (Pn,n) $n,mSa,b- 



In this case, we can solve analytically for the absorption 
because this restriction leads, from Eq. (|62|) to 



(0)o,6,c,d / q 



(83) 



— Xni,n 2 ',ri2,ni (^i ^n) 5 a ,d5b,c5n 1 ,ni5n 2 ,n s - 

Now, at q = the only nonzero Fock interactions in the 
GRPA equations are given in Eqs. (|49M51|) while the only 
Hartree interactions that need to be considered are those 
of the form 



HlloXo (0) , flf^f,i (0) , H^l, (0) , //^Ao (0) • (84) 

These intcrlayer Hartree interaction involves combina- 
tions of the form e~ qd /q that give a finite contribution at 
q = (and also a diverging contribution that is cancelled 
by the other terms) . It follows that Eq. (loTfl) gives for the 
GRPA response functions Xi'o'i'o ( w ) = Xo'i'o'i ( w ) = 
and 



a.a.c.c i \ 

Xo 1,1 o M 



(0)a,a,a,a 

Xo, l.i.o 



(w) 5 a ,, 



1 + n A i, i,o,o ( u )X ,i,i,o M 



(85) 



i.e. valley and spin must be conserved in an optically 
active electronic transition. Since 



(0)a,a,a,a , % _ 

Xo, 1,1,0 l w J — 



Wo) - (Pl.l) 



u + l 5-(A a + uZ: a -u? : n /h 



(and a similar expression with 
we have easily 



(0)a,a,a,a 



a, a, a, a / \ 

Xo i i,o M = 



1 for v w ' ' 
1 101 Xi, 0,0,1 

/ a,a\ / a,a\ 

Wo) - (Pl,l) 

Ul ■' 



and 



a. a. a. a / \ 

Xi o o,i M = 



A2 + ^ </>;;?>] /« 

/ a,a\ / a,a\ 

W'l) ~ W,o) 

.. + W+[A0 + ^(p?;?>] /ft 



(86) 

M), 

(87) 



tii r a, a, a, a / \ j a, a, a, a / \ . 1 

Ihe functions X 1 1 o ( w ) ana - Xi o o l ( w ) are the re ~ 
sponse to the two circular polarizations of light. The 
absorption in an incoherent phase is finally given by 



El 



E 



<Po°;o°>-<Pm» ( 89 ) 



(K) 2 5 



fa - A° - ^ <p£) 



with a — x,y and A° given by Eq. (1761) . In the numerical 
calculation, we introduce a small Landau level width in 
order to get a finite value for the optical absorption. 

We can follow the same type of reasoning to show that, 
in a phase with no orbital coherence but with possibly 
layer and/or spin coherence, the functions Xi'o'i'o ( w ) = 
Xo'i'o'i i 1 ^) = ano - the absorption depends again only 

a.a.b,b / \ j a.a.b,b / \ T ,1 • • i , , 

on Xo 1 1 o ( w ) ano - Xi o o i ( w ) ■ ^ n this special case, the 

a,b,c,d 



(82) equation of motion for Xo'l'i'o ( w ) ^ s 



a,a,b,b / \ (0)a,a,b,b , \ 

Xo,i,i,oH = Xo,i,i,o M 



(90) 



2^ ^0,1,1,0 M A l,l,0,oX ,l,l,0 M 



and a similar expression with 0^1 for Xi'o'o'i ( w ) 



IV. PHASE DIAGRAM OF THE C2DEG 

At zero bias, the QHF states follow a set of Hund's 
rules: the spin polarization is maximized first, then the 
layer polarization is maximized to the greatest extent 
possible and finally the orbital polarization is maximized 
to the extent allowed by the first two rules^. In this 
section, we study the phase transitions that occur when 
a finite bias (or transverse electric field) is turned on. 



A. Types of phases 

Fig. [5] shows our numerical result for the phase dia- 
gram of the C2DEG as a function of an applied transverse 
electric field E = Ab /ed for B = 10 T and k = 5 and 
for all integer filling factors v £ [— 3,3] . We indicate the 
eight non-interacting levels by horizontal lines and num- 
ber them according to the scheme indicated in the top in- 
set. Note that the lines are only offset vertically for clar- 
ity. Their position does not reflect the true ordering of 
the energy levels which changes with bias. We name the 
phases /„, L„, O v and SL V according to the type of coher- 
ence that is present: incoherent, layer-coherent, orbital- 
coherent or spin-layer coherent respectively. When there 
is more than one incoherent phase at a given filling fac- 
tor, we use the notation I* for the second phase, /** for 
the third phase and so on. The critical electric field for 
the transition between two phases is indicated by to 
i?( 18 ) and is in units of mV/nm. A circle on an energy 
level represents a fully filled level while an ellipse that 
connects two levels indicates a coherent superposition of 
these two states. We list in Fig. [5]some properties of each 
phase: spin polarization S z , number of Goldstone modes 
(G), number of collective modes gapped at the Zeeman 
energy (Z) and number of peaks in the optical absorption 
spectrum (A). The transition between a coherent and an 
incoherent phase is continuous while a transition between 
two incoherent phases is discontinuous. 
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L ,:S =1/2 (1G,2Z,2A) 



E =0.11 
L 2 :S 2 =1 (1G,2Z,2A) 

E (2) =0.35 

L,:S =3/2 (1G,2Z,2A) 

• • — ■,— 

E""=0.23 

VS =2 (0G,2Z,0A) 

• • 

L,:S 2 =3/2 (1G,2Z,2A) 



E'"=0.11 
L 2 :S =1 (1G,2Z,2A) 

E IU, =0.35 
L,:S =1/2 (1G,2Z,2A) 



I ,:S =1/2 (0G,1Z,1A) 



I,:S =1 (0G,1Z,0A) 



K + ,+,l 


K.,+,1 


K + ,-,l 


K,-,l 


7 


3 


8 


4 


5 


1 


6 


2 


K + ,+,0 


K,+,0 


K + ,-,0 


K,-,0 



E =459 



I 2 :S Z =0 (0G,0Z,1A) 



I,:S =3/2 (0G,2Z,1A) SL,:S a e [3/2,1/2] 
1 (1G,2Z,2A) 



E =11.94 



E'=12.! 



E"'=12.1S 



SL:Se[2,0] 
(1G,2Z,2A) 



E"=0.23 



E I8I =13.62 

I i: S 2 =3/2 (0G,2Z,1A) SL :S g [3/24/2] 
• (1G,2Z,2A) 



E ll0 =14.44 
L,:S 2 =1 (0G,1Z,0A) 
• • • 

L:S =1/2 (0G,1Z,1A) 



E""=14.85 



I,:S =1/2 (0G,1Z,1A) 



I^:S 2 =0 (0G,0Z,0A) 



I^S =1/2 (0G,1Z,1A) 



O i: S 2 =l/2 
(1G,1Z,0A) 



E" 2, =135 



+nonuniform states 
> 

E II3I =500 



E' ls '=1097 

0,:S =1/2 
(1G,1Z,0A) 
+nonuniform states 



E II7I =135 



E" =500 



£ :S=l/2 
(0G,1Z,1A) 



I,:S 2 =0 (0G,0Z,1A) 



(0G,1Z,1A) 



E (mV/nm) 



FIG. 5: (Color online) Phase diagram of the C2DEG in Landau level N = at B = 10 T and for function of 

the transverse electric field beween the layers for integer filling factors v G [—3,3] (from top to bottom). A filled (red) circle 
represents a filled state while a filled (blue) ellipse indicates a coherent superposition of two levels. The numbering of the levels 
is indicated in the inset at the top right of the figure. The E^'s indicate the critical perpendicular electric field (in mV/nm) 
required for the transition between two phases. Also indicated for each phase are the spin polarization S z , the number of 
Goldstone mode (G), of modes gapped at the Zeeman energy (Z) and the number of peaks in the optical absorption spectrum 
(A). The mention "+ nonuniform states" for the Oi and O3 phases signals that this portion of the phase diagram is further 
subdivided into uniform and nonuniform states as indicated in Eq. (|109p . 



Incoherent phases I v 



The l v phases have no coherence of any kind and so 
L||,5||,0|| = (the parallel component is in the plane 
of the bilayer). Each level is either full or empty and 
so L Z ,S Z and O z vary from one phase to another. The 
state corresponding to a specific diagram is easily read 
from Fig. 5. For I2, we have for example 



n 

X 



C 7,X C 5,X 



|0>, 



(91) 



and the order parameters (ps^) — {pij) = 1. We include 
the phases I± 2) I** and 1% in the phase diagram only to 
make it more complete. Indeed, the bias can't produce 
any more transition after these states. But, the critical 
bias needed to reach these states is well outside the limits 
of validity of our two-band model. 

Using Eq. P5|) . the Hartree-Fock energy of two adja- 



cent incoherent phases are readily compared to extract 
the critical biases. We find, with = edE^ , 



1+2 -> II, : A 



(15) 





A z + Ci + | A c 



+2 • 



(92) 
(93) 



If we ignore the SL V phases at v = 0,±1, we find that 
the transition between the two incoherent phases occurs 
in the middle of the SL V phase (see Sec. IV. A. 4 below) 
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i.e. at 



where the critical bias is, to order (d/£) , 



i_i -> ilj : A_x = A z + 2-a, 

_ A z + 2|a 
1-/3 ' 
A z + 2§a 



In -► In* : A 



Ji -> If : A 



1 - 2/3 



(94) 
(95) 
(96) 



Note that these last three results are independent of 
Coulomb exchange corrections and so of screening cor- 
rections. (The capacitive term comes from the Hartree 
self-energy and is not screened.) 

From Eq. ([55]). we find a critical electric field E c = 
Ao/ed 4. 75 [T] mV/nm when k = 1 for the transi- 
tion Iq —> Jq. This critical field depends linearly on the 
magnetic field, in agreement with the experiments 12 . Ex- 
perimentally, however, the slope isi£ 11 mV nm _1 T _1 
ori£ 12.7 mV nm" 1 T" 1 or 14 12 - 18 mV nm" 1 T" 1 and 
thus larger than the HFA value. Moreover, experiments 
measure an offset of E c ss 20 mV/nm at B = T— . 
This offset can't be captured by our HFA which is only 
valid at sufficiently large magnetic field where Landau 
level mixing can be neglected. Apart from the extra j3 
corrections, Eqs. (I94H96[) are identical to those given by 
Gorbar et al*i£. 



2. Layer-coherent phases L„ 

The second type of phase, L v , has layer coherence be- 
tween two states with the same spin and orbital indices 
and so i|| ^ but S\\ , 0\\ =0. The tilt angle of the pseu- 
dospin vector L varies with bias in this phase but O z 
and S z are constant. An example is phase L_ 3 which is 
described by 



i*£-a>=n( ac U +6c U) 

X 



(97) 



where the coefficients a and b depend on the bias and are 
related by \af + \b\ z = 1. With increasing bias, a — > 1 
and b — > continuously. The level populations and the 
coherence in £_3 are given by 



(pB,s) = \a\ 2 , (pi.i) = \b\ 2 , 
(pi,s) = (P5,i>* = ab*. 



(98) 
(99) 



Fig. [6] shows how these variables depend on the trans- 
verse electric field for the similar phase L-\. The popu- 
lations of the coherent levels vary linearly with the bias 
in all L v phases with the exception of L±2 where the 
variation is not exactly linear. In for example, 



(P5,5 




(100) 



3(1) 



-3(1) 



A« = Af 



i-i (3 ) -> /- 1( 3) : AW = A^ 6) 
7 fir f<P" 



1(7) 

(102) 



These critical biases all scale with the magnetic field as 

£3/2 

In phases L±2, there is a layer coherence in orbitals 
n = and n = 1. Phase L_2, for example, is described 
by the state 

|*l_ 2 ) = U{ ac lx+^{ x ) (103) 

X 

x («'4,x+Y4j io>. 

The critical bias A^ 2 ' = A^ 14 ) for the transitions L±2 — > 
I±2 has a complicated analytical expression that we do 
not reproduce here but the numerical values of the critical 
electric field for B = 10 T and k — 5 is indicated in Fig. 

El 

Phases with layer coherence occur in a small range of 
bias and at very small bias because the interlayer sepa- 
ration d/£ — 0.013"\/B is very small in bilayer graphene 
and so is the capacitive energy. In semiconductor bilay- 
ers, d/l can be of order unity and interlayer coherence 
can survive to a much higher bias^i. 



■s 0.4 



Phase L, 




" < P7.7> 








<P 3 ,3> 







0.1 0.2 
Electric field (mv/nm) 



FIG. 6: (Color online) Variation of the populations and in- 
terlayer coherence with the transverse electric field in phase 
L_i. 



3. Orbital-coherent phases O v 



For A B > A^ 2) 



f 17) 

= A^ , the ordering of the energy 
levels n = 0, 1 is reversed (this change is not shown 
in Fig. [5J. When this happens, the kinetic energy is 
minimized by filling level n = 1 before n = 0. How- 
ever, this increases the Coulomb exchange energy because 
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x ti + i.i (°) < X oao.o (°) ■ The C2DEG optimizes its en- 
ergy by creating a coherent superposition of n = and 
n = 1 with the same valley and spin indices. We use the 
notation O v for such a phase. An example is phase O3 
which is described by 

l*o s ) = Y[(c4,x + telx)4,x4,x ( 104 ) 

X 

XC 7,X°5,X C 3,X C 1,X |0) • 

In this state, L|| = Sn = and L z and are constant. 
It is now the tilt angle of the pseudospin vector O that 
varies with bias. Orbital coherence begins when the bare 
energy of the state |if_,±,0) is equal to that of state 
\K-, ±, 1} at v = 1, 3. This occurs when 

I* -> Ox : Ag 2) = |, (105) 

7 3 -► 3 :Ag 7) =Ag 2) , (106) 

and the critical bias does not depend on the magnetic 
field, Coulomb interaction or on the value of the dielectric 
constant. We find Ci/fi = 46 meV i.e. E — 135 mV/nm 
which is in the range of validity of the two-band model. 

The orbital phase survives until a transition to an in- 
coherent phase occurs at the critical bias 

1 -+ I*,* : a£ 3) = j (Ci + |Ac) , (107) 
3 -> I 3 * : a£ 8) = a£ 8) . (108) 

Since is the onset of the orbital phase, we see that 
the range of existence of the orbital phase scales as 1/ y/B. 
Part of this range is outside the limit of validity of our 
model. 

Fig. [7] shows how the populations and coherence vary 
with the transverse electric field in phase O3. The same 
behavior is found in phase 0\. We remark that in pre- 
vious work where the spin degree of freedom is frozen^, 
the orbital coherent phase occurs at v = — 1 and v = 3. 




Electric field (mV/nm) 

FIG. 7: (Color online) Variation of the populations interlayer 
coherence with the transverse electric field in phase O3. 



The orbital phase exists in a large range of bias. As 
we explained in Sec. III.E, a finite orbital coherence im- 
plies a finite density of electric dipoles in the plane of the 
layers. The orientation of these dipoles can be controlled 
by an electric field in the plane of the layers^. On of 
us has studied in detail the interesting properties of this 
stat o 25 ' 26 . For example, the collective mode associated 
with the orbital coherence is highly anisotropic. This 
mode softens at a finite wave vector in the direction per- 
pendicular to the dipoles when the bias is increased. This 
suggests a transition to a charge-density-wave state. In 
the Hartree-Fock approximation, it was found that this 
transition is preempted by a transition to a crystal phase 
with one electron per site and a Skyrmion-like pseudospin 
texture of the orbital pseudospin at each crystal site. As 
the bias is increased, the crystal state is followed by a 
helical state where the orbital pseudospin rotates along 
one spatial direction. In both phases, the total electronic 
density is modulated spatially but the local filling fac- 
tor is not. By further increasing the bias, the crystal 
state and then the uniform states are recovered. The 
critical electric fields for the transition to the uniform 
(UP), Skyrmion crystal (SKP) and helical phases (HP) 
are given b y 26 i 32 : 



134 < A B < 145 mV/nm 


UP 


145 < A B < 189 mV/nm 


SKP 


189 < A B < 450 mV/nm 


HP 


450 < A B < 494 mV/nm 


SKP 


494 < A B < 505 mV/nm 


UP 



The phase diagram in 0\ and O3 is symmetric with re- 
spect to the center of the helical phase. 

We remark that this sequence of phase transitions is 
similar to that observed in a thin film of the helical 
magnet Feo.5Coo.5Si when a perpendicular magnetic field 
is increased^ 3 -. It has been shown 2 ^ that the Hamil- 
tonian of the C2DEG in the orbital phase contains a 
Dzyaloshinskii-Moriya (DM) interactional that is respon- 
sible for the rotation of the pseudospins. Its origin in 
bilayer graphene is purely Coulombic while the DM in- 
teraction comes from spin-orbit coupling in Feo.5Coo.5Si. 

We do not find any sign of instability in the collec- 
tive mode dispersions for the other phases in Fig. [5] 
However, we remark that phases with lower energy than 
those considered in this figure are possible. In order to 
establish the phase diagram of the C2DEG, we choose a 
set of possible ground states and compare their energies. 
This does not ensure however that the true ground state 
is amongst the states that we have chosen to compare! 
For the incoherent states, this is not a problem because 
there is a finite number of states to compare. But for the 
nonuniform states, the number of possible ground state 
is enormous. 
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4- Spin-layer- coherent phases SL V 

The fourth type of phase has coherence between two 
states with the same orbital index but different spin and 
layer indices. We use for these phases the notation SL V . 
An example is phase SX_i where the ground state is 

1**0 = II { ac lx + bc{ x ) 4 x c{ x |0) . (110) 

X 

Because the coherence is now between two states with 
different spin and layer indices, i.e. \K +1 — , 0) and 
\K-,+,0) in we cannot describe the change 

with bias as the tilting of one of the pseudospin L, S or 
O. In fact this state has Ln, Sn = On = 0. Both L z and 
S z vary with bias however. This phase is characterized 
by the order parameter (p6,i) • 

The variation of (p6,e) i (pi.i) an d (p6,i) with bias in 
this phase is similar to that shown in Fig. [7J The critical 
biases for the beginning (A^ (n)) and end (A e (n)) of the 
SL±i phases are given by 



(111) 
(112) 

(113) 
(114) 

(115) 
(116) 




with the exception of the phases SL V where it varies 
continuously between the two numbers indicated. The 
biggest change in S z and L z occurs at filling factor v = 
where the C2DEG makes a transition from a fully spin 
polarized and layer unpolarized gas (S z — 2fi, L z = 0) 
at small bias to a spin unpolarized and layer polarized 
gas (S z — 0, L z = 2) at large bias. For v = ±1, the 
SL V phase interpolates between S z = 3fi/2,L 2 = 1/2 
and S x = h/2 : L z = 3/2. The only jumps in S z occur at 
the transitions I ±2 — > I± 2 where the system goes from a 
spin polarized to a spin unpolarized state. 



C. Transport gaps 

Another quantity that is accessible experimentally is 
the transport gap A which is defined by the difference in 
energy between the first empty state and the last filled 
state of the Hartree-Fock Hamiltonian. It was shown 
previously^ that the gap at zero bias follow the hierar- 
chy A„ = o > A„ = ±2 > A„ = ±i.±3. This implies that the 
first plateau to appear when the magnetic field is turned 
on has <7 xy = 0. At larger field, the a xy = ±2 plateaus ap- 
pear and at still larger field, the a xy — ±1, ±3 plateaus. 
This is indeed what is seen experimentally^^. 

Fig. [5] shows the Hartree-Fock gaps as a function of 
the transverse electric field for the different phases of the 
C2DEG. For this figure, we have taken B = 10 T and k = 
5. In some phases, one or more level crossing occurs that 
change the behavior of the gap. This is clearly visible for 
I ±2 in Fig. M (a) and for J_i in Fig. [8] (c). 

With the exceptions of the phases where coherence oc- 
curs in two levels (L±2, SLq) 7 it is possible to obtain a 
simple analytical expression for the gap. We list these ex- 
pressions below. When one or more level crossings occur, 
we use the notation I„ ,lj> S \ ... to denote the differ- 
ent behaviors of the gap and A^ f or values at 
which the level crossings occur. 

In the incoherent phases with v = ±1, ±3 : 



These critical biases scale linearly with the magnetic 
field. The phase SL V is the ground state in a small range 
of bias of the order (d/£) 2 a which is approximately 0.1 
meV for k = 5. 

A sufficiently large bias is necessary for spin down 
states to cross the spin up states and produces a SL U 
phase. The SLq is special because it involves coherence 
in both n = and n — 1. As for L±2^ the exact critical 
bias in this case has a complicated analytical expression 
which we do not reproduce here. 



j t r*(l) 



f3A B + Ci + -A c , 



/ 3 *,/r : A = /3A B - Ci + -A c . 



/_i,/r,7 3 : A = -(3A B + Ci + -A c , 



'?> : A*_, =A Z + |A C , 



(117) 
(118) 

(119) 
(120) 



B. Spin polarization 

We indicate for each phase in Fig. [5] the spin polar- 
ization S z . The polarization is constant in all phases 



where, for v = — 1 : 

a(i)-(2) _ A 2 - Ci + |Ag 



B 



P 



(121) 
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For the incoherent phases with v = — 2 



by 



r( 2 ) 



A = (l-/3)A B -Ci-2^a+^A c , (122) 
A = -pA B - Ci + A z - 2^a + ^A c , (123) 



r 2 : A = /?A S +Ci-A z + -A c , 



(124) 



with 



A 



U)-(2) _ 
B ~ 



Az + 2-a, 



(125) 



while for v = 2, 



4 1 ' = A = (l-/3)A B -Ci-2^a+^A c , (126) 
: A = /3A B +A z -Ci + ^A C , (127) 



r(3) 



A = Ay + -A r 



(128) 



: A = -/?A B + A z + Ci + ^ Ac, (129) 



with 



A (i)-( 2 ) _ 2 ? a + A, 
5 " 1-2/3 : 



A 



(2)-(3) _ Cl 



A (3)-(4) _ Cl + jA C 
/3 



(130) 
(131) 
(132) 




(137) 



(138) 



where A±i are defined in Eqs. (|94I96[) . The gaps are 
twice as big in phases Iq, I ±2 than in phases I±\, I±$ and 
vary more rapidly with bias in the former than in the 
latter. The presence of the SL U phase smoothens the 
jump of the gap in the transition from I±i to 1^. 

The gap is independent of the bias in the orbital phases 
Oi and O3: 



0i,03 : A 



(139) 



For the L v phases, the gaps at zero bias are given ap- 
proximately by 



A c , (140) 




A « -Ci - 2 ( - ] a 



(142) 



For v = 0, we find: 



la ■ A = - (1 - 0) A B - Ci + A z + -A c , (133) 



Io (1) ■ A = (l-/?)A B -Ci-Az 

A 5 A 

_ 4 -« +i A c , 
I*o (2) : A = (1 - 2/3) A B - A z 



with 



A 5 A 



A 



(l)-(2) = Cl 



(134) 
(135) 



(136) 



The gap changes rapidly in the SL±\ phase while it is 
almost independent of the bias in phase SLq. Its value 
in the middle of the SL±i phases is given approximately 



The correction (d/l) is very small and the gaps at zero 
bias for L±\ and L±3 are almost equal. (The difference 
comes from the fact that the coherence is not in the 
same orbital in L_3,Li and L_i,i3.) These gaps are 
not shown in Fig. [5] because the corresponding phases 
occur at very small biases. The gap increases (decreases) 
with bias in phases L_3,Li (L-\,Ls). It is almost con- 
stant in L±2- The main contribution to all gaps is the 
Coulomb exchange interaction. 

With Ci = 0, our gaps for phases I2,Io,Iq agree with 
those of Gorbar et alJ^ if screening is neglected in their 
calculation. For ii,Jx and -Z3, however, our exchange 
correction is Ac/2 which is consistent with Ref. Il3l while 
Gorbar et al. have 3 Ac/8. 

At A B — 0, we find for Iq the gap A = 62 meV for 
B = 2 T and K — 2 while, with static screening, Gorbar 
et al. find ps 5 meV. Similarly, for phase I\, the HFA 
gives A = 26 meV while the result with static screen- 
ing is ss 2.5 meV. Static screening leads to a reduction 
of the gap by a factor of at least 10. Dynamical screen- 
ing and Landau mixing corrections, however, increase the 
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gaps calculated with static screening by a factor of two 
to thresh. As for the behavior of the gap with bias. 
Fig. 5 of Gorbar et al^ shows that, with screening, the 
gap of the phase ^3 increases with the electric field even 
when the correction /3As is neglected. The slope is ap- 
proximately 0.1 nm-C for B = 2 T and k — 2. The 
(unscreened) HFA predicts a slope of j3d = 0.006 nm-C. 
When both screening and /3As are considered, the gaps 
for I±i, I±3, I±i will probably increase with bias (con- 
trary to the behavior illustrated in our Fig. [5£b)) but 
the rapid change of the gap in the SL±\ will still be 
present. 

The energy gaps obtained from local compressibility 
measurements on suspended bilayer graphene by Martin 
et al& are of size A u=0 w 1.7B [T] meV, A„ =±2 1.2B 
[T] meV and A v= ±i « 0.LB [T] meV (with less data 
points in this case). A more recent transport experiment 
by Velasco et al 10 on suspended bilayer graphene with 
a higher mobility reports a larger gap A„ = o w 5.5B [T] 
meV. The measured gaps scale linearly with the magnetic 
field contrary to the HFA prediction. In fact, Gorbar 
et aZ.fi^ have shown that a linear scaling is obtained if 
static screening is considered. (One set of experiments 
at higher magnetic field reported gaps that scaled as \f~B 
however—.) 

Experiments^ show that O 'xy CGclSGS to be quantized 
at v = 0, 1 in the region corresponding to the SL V phase 
and at v — 2 and v = 3 in the region around zero bias. 
A possible explanation is that the conductance quantiza- 
tion is broken by disorder in the regions corresponding 
to a minimum of the gap^£. However, that argument 
does not seem to work at v = 1 where the screened HFA 
gap increases continuously with the bias and is also not 
compatible with our unscreened result. 
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D. Collective modes and optical absorption 

Each phase of the C2DEG is characterized by a set of 
collective excitations. The number of dispersive modes 
when to levels are filled is to (8 — to). We have calculated 
the dispersion relation of these modes using the GRPA 
described in Sec. IIII. In the limit q — > oo, the vertex cor- 
rections vanish and the response function x ~ ► X° where 
X° is the response function evaluated in the HFA. Thus, 
in this limit, the collective mode frequencies must corre- 
spond to transitions between a filled and an empty eigen- 
state of the Hartree-Fock Hamiltonian. At finite value of 
q, some modes mix together and it becomes difficult to 
identify their character (layer, orbital, spin transitions, 
etc.). Our numerical results are summarized in Fig. [S] 
for the coherent phases and in Fig. [TU]for the incoherent 
phases. 



FIG. 8: (Color online) Variation of the transport gap with 
the perpendicular electric field in the different phases of the 
C2DEG in Landau level N — 0. The kink in the behavior of 
I ±2 are due to level crossings. A value of n = 5 is assumed 
for the dielectric constant and B — 10 T. 



1. Goldstone modes 

The coherent phases sustain one gapless (Goldstone) 
mode. The number of Goldstone modes is indicated for 
each phase in Fig. [5] For example, in phase £_3, this 
mode is due to the fact that the layer pseudospin L can 
rotate freely around the z axis. The same situation oc- 
curs for the coherent phase 0„ where again the orbital 
pseudospin O can rotate freely around the z axis. Phases 
L±2 support coherence in both n = and n — 1 orbitals 
and we can define a layer pseudospin Lq for n = and Li 
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FIG. 9: (Color online) Dispersion of the collective modes in the coherent phases. The (blue) arrows point to the modes which 
are active in optical absorption experiments. All modes are evaluated at B — 10 T with k = 5. The bias in units of € 2 /k£ = 35.6 
meV (corresponding to a frequency v — e 2 /hn£ = 8.6 x 10 12 Hz) is respectively: As = 0.0005 for L_3,Li; Ab = 0.001 for 
L 3 ,L-i;A B = 0.002 for L_ 2 ; A B = 1.26 for 3 ; A B = 1.115 for SX_i;A s = 0.14 for SLi and A B = 0.126 for SL . In SL 
the middle line in each group of three dispersive curves contains two modes which are very close in energy. 



for n — 1. The Goldstone mode in this case correspond to 
an in-phase rotation of both pseudospins. Alternatively, 
we can see this mode as a fluctuation of the relative phase 
of the two order parameters (pifi) and (^03,7) in L_2 or 
(P2,e) and (p 4)8 ) in L 2 . 

The dispersion off all gapless modes (with the excep- 
tion of the gapless mode in the orbital phase) is lin- 
ear in wave vector at very small wave vector i.e. for 
1^ ~ d/L Phases L± 2 have the same collective mode spec- 
trum. The dispersion of the Goldstone mode in phases 
where the coherence occurs in the orbital n = 1 (i.e. 
L—i, L3, SLi) has a roton minimum while there is none if 
the coherence occurs in n = (i.e. L_3, L\, SL—±). This 



is due to the particular form factor for n — 1 involved in 
the Coulomb matrix elements (see Eqs. (|44H47[) ). Phases 
L—2i L 2 , SLq contain coherence in both n = and n = 1 
and a small shoulder appears in the dispersion. 

For phase L v , the Goldstone mode is the famous layer- 
pseudospin-wave mode which has been extensively stud- 
ied in semiconductor bilayer— at filling factor v = 1 and 
detected experimentally^. In semiconductor bilayer, this 
mode becomes soft at a finite wave vector as the separa- 
tion between the layers is increased (around djl s» 1). In 
bilayer graphene djl « 1 and the layer-coherent phases 
are stable. The only instability in the collective modes is 
seen in the orbital phases 0\ and O3. 



17 



For (i = 0, the Goldstone mode of 7_ 3 (and L{) has a 
quadratic dispersion at zero bias and becomes unstable 24 
at finite bias. (We have checked that these conclusions 
remain valid if £1 is finite but small.) A consequence of 
this instability is that the C2DEG is expected to go from 
a smectic (non-homogeneous) phase at low temperature 
to an isotropic phase at higher temperature. The smectic 
phase would lead to anisotropic electrical transport. 

The dispersion of the gapless orbital pseudospin-wave 
mode was studied in detail for the uniform phased as well 
as for the crystal and helical phases 2 ^ In the uniform 
phase, it has a strongly anisotropic dispersion: linear in 
the direction of the orbital pseudospins and quadratic in 
the other directions i.e. 



w(q) - V2(/3A B -Ci)^|sin(0 q )|, 



(143) 



w(q) = ^v^O8A B -CiM,(0q = O,7r), (144) 

where f? q is the angle between the wave vector and the x 
axis. The Goldstone mode softens at a finite wave vector 
qi 2 in the direction perpendicular to the orbital pseu- 
dospins at a bias = 58.8 meV (E = 173 mV/nm). 
This suggests a transition to a charge-density- wave state. 
As we explained above, this transition is preempted by a 
transition to a crystal phase at = 47.7 meV. 



2. Spin-wave modes 

The number of spin-wave modes gapped at Az is in- 
dicated for each phase in Fig. [5] 

All coherent phases with the exception of 0\ , O3 have 
two modes gapped at the Zeeman energy Az at q = 0. 
Their degeneracy is lifted at finite wave vector. Because 
of the occupation of the levels in 0\ and O3, only one 
intralayer spin-flip transition is possible in theses phases. 

The incoherent phases can have 0, 1 or 2 modes gapped 
at Az- To be gapped at Az, these modes must involve 
transitions within the same valley and orbitals. From 
Fig. [5j it is easy to see that no intralayer and intraor- 
bital spin-flip transition is possible for 7q and I± 2 and 
that the occupation of the levels permit only one such 
mode in i±3,J±i- In I±2, two transitions seem possible 
but they are degenerate and the coupling between them 
leaves one mode gapped at Az and the second mode has 
its frequency renormalized. The same mechanism oper- 
ates in phases I±i,Iq resulting in two modes gapped at 
A z . 



3. Optical absorption 

In the absence of Coulomb interaction, the dynamical 
conductivity has intra-octet peaks at the bare gap en- 
ergy = - E^ Q = Ci + S,/3A B in addition to the 
inter-Landau-level peaks which do not appear in our cal- 
culation. Fig. Illl shows the absorption in different phases 



when Coulomb interaction is considered. The number of 
absorption peaks is also indicated for each phase in Fig. 
[5] and, in Figs. 191101 we identify by (blue) arrows the 
modes that lead to optical absorption. 

In the incoherent phases, optical absorption is possible 
only for transitions that occur between states with the 
same valley and spin indices but different orbital indices. 
The incoherent phases have one or zero absorption peak. 
The latter case occurs when both levels n = 0, 1 with the 
same layer and spin indices are filled (phases I ±2, Iq, Iq ). 

When level n = is filled, the absorption, as 
shown in Sec. IILg is exactly at u> = (Ci +£,/3Ab) /h. 
When level n = 1 is filled, the absorption is at u> = 
(Ci — (3Ab + Ac/4) jh i.e. affected by exchange correc- 
tions. The former case applies to most of the incoherent 
phases in our phase diagram. The latter case applies to 
phases 7-j** ,7^ an d 7g which occur at very high bias and 
are outside the limits of validity of our two-band model. 
In those phases, level n — 1 is filled because it is be- 
low n = in energy. In phases I ±2, the two allowed 
transitions have the same energy and the intensity of the 
absorption peak is doubled. 

Fig. HU (a) shows the absorption as a function of fre- 
quency in phase 7_3 at finite bias. The absorption is 
concentrated in one strong peak at the frequency 



v = (C 1 +^A B )/h 
= (9.8 + 7.6£A B 



(145) 



k(. 



x 10 11 Hz 



at B = 10 T and k = 5. The absorption frequency varies 
widely with bias. For example, at the onset of the transi- 
tion from 7* to Oi, the frequency v — > (see Eq. (|105|) ) 
while j/«9.8xlO n Hzat the onset of the I± phase. By 
contrast, if the absorption occurs in layer 7C|_, the mini- 
mum frequency is v = 9.8 x 10 11 Hz since the absorption 
frequency must increase with bias in this case. 

In the Oi,C?3 phases, the Goldstone mode has orbital 
character (i.e. electric dipole fluctuations) but does not 
lead to absorption at finite frequency. 

Interestingly, the coherent phases L v and SL V show 
two absorption peaks (see Fig. QTJ). The second peak 
in 7±i,7±3 is extremely weak however and is absent 
at zero bias. The two weak peaks in L±2 are also ex- 
tremely weak and disappear at zero bias. In 7_3,7i, 
the first absorption peak at zero bias is exactly at 
huj = (1 while for 7_i,L3, the frequency is at fiw = 

Ci + {Kfito,o - *oao,o " KCiA + I 2 which is 

slightly shifted from Hut — £1. 

The layer eigenstates of the Hartree-Fock Hamiltonian 
at zero bias are the symmetric (S) and antisymmetric 
(AS) combinations of K + and 7T_. Using Eqs. ([90]) 
and ([Hll), we can show from the GRPA equations that 
the response functions that enter in Eq. (f5Tj) for the 
absorption in the absence of orbital coherence and at 
zero bias are of the form ^^(^-iM^U^-i) (w) 
or X ^UH^)M>^)M^)[^ where a? now stands 
for the S or ^45 layer combinations. That is, the layer 
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FIG. 10: (Color online) Dispersion of the collective modes in the incoherent phases for B = 10 T and k = 5. The bias is in 
units of e 2 / k£ = 35.6 meV corresponding to a frequency v = e 2 /hud. = 8.6 x 10 12 Hz. The (blue) arrows point to the modes 
which are active in optical absorption experiments. The bias in units of e 2 ' jnl is respectively: As = 0.01 for Iq, i±i, I±2, 1 ±3 
and As = 1.0 for Iq,I±i. For 7±3,/±2, i±i only the first four, three, and six modes respectively are shown. For 7o, each of the 
four branches contains four modes which are close in energy. In Iq the middle line in each group of three curves contains two 
modes. 



combination is conserved in the absorption at zero bias 
and in the absence of orbital coherence. It follows 
that, in phase £-3, the state \S,+, 0) is filled and only 
the transition \S, +,0) — > |5, +, 1) is optically active. 
For L_2, no transition conserving the valley index (S 
or AS) is possible and for L_i the only allowed tran- 
sition is between \AS, +,0) — > \AS,+,1) since levels 



\S, +, 0) , \S, +, 1) , \AS, +, 0) are filled. The same argu- 
ment applies to the spin down states. 

At finite but small bias, the S and AS layer combina- 
tions are replaced by bonding and anti-bonding combi- 
nations i.e. \B,+,0) = a\K_,+,Q) +b\K+,+,0) and 
\AB,+,Q) = -b\K-,+,0) + a\K+,+,0) for example 
where a and b depend on the bias and on the orbital 
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index n. The strong peak in the absorption corresponds 
to the transition \B,+,0) — > \B, +, 1) for L_ 3 and to 
|AB,+,0) \AB,+,1) for L_i. By contrast to the 
zero bias case, the absorption given by Eq. (fT9"|) con- 
tains response functions which are not just of the form 
x (uj,^)M^)M^Ua ( ,^) (w) with = B AB SQ that 

transitions that do not conserve the layer combination 
B,AB are weakly optically active 3 -. A second peak ap- 
pears in the absorption spectrum which corresponds to 
the transition \B,+,0) — > \AB,+,1) for L-3 and L-%. 
The two weak peaks in L_2 come from the transitions 
\B, + , 0) -> \AB, + , 1) and \B, +, 1) -> |AB, +, 0) . 

Fig. [TT] (b) shows the absorption in phases SXq an d 
SL—i. In SX_i, we find by analyzing the eigenvectors of 
the modes involved in the optical absorption that the 
strong peak corresponds to the transition \B, r\, 0) — > 
l-B, 77, 1) and the weak peak to the transition \B, r\, 0) — > 
|AB,77, 1) where |-B,ry,0) = o|if_,+,0)+6|if+, -,0) and 
j AS, 77,0) = — b |iir_, +, 0) +o|Jf+, -,0) and a, b depend 
on the bias and on the orbital index n. The layer com- 
bination here are between two states with opposite spin 
orientations. 

The electromagnetic absorption in the crystal and he- 
lical phases are much more complex and was discussed 
previously 26 . In the helical phase, for example, the ab- 
sorption depends on the orientation of the polarization 
of the electromagnetic wave in the x — y plane. 

The collective mode dispersions computed here for 
the intra-Landau level [N = 0) transitions should be 
combined with the dispersion of the inter-Landau level 
magnetoexcitons computed by Sari and Toke 2 ^ and by 
Shizuya^ to get a complete picture of the absorption for 
the C2DEG in bilayer graphene. 



V. EFFECT OF AN IN-PLANE ELECTRIC 
FIELD 

In this section, we consider the effect of adding a uni- 
form in-plane electric field to the Hartree-Fock Hamilto- 
nian of Eq. (|4Tj) . The coupling between the C2DEG and 
this external electric field is given by 



H E = -d-E h (146) 

where d is the total dipole moment of the electron gas 
given by Eq. (|70|) . 

The main effect of E|| is to induce orbital coherence. 
We consider here the case of v = —1 but a similar ef- 
fect occur at other filling factors and will be discussed 
elsewhere^. We set Em = —E x. In the ground state, 
the electric dipoles are aligned with E||. 

With finite Em, the SL_i state is replaced by a state 
with spin, orbital, and layer coherences i.e. SOL—\. This 
is represented by the inset in the top-left corner of Fig. 
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FIG. 11: (Color online) Electromagnetic absorption for B — 
10 T and ft = 5 in phases (a) L_ 3 (A s = 0.018 meV), Is 
(A B = 0.36 meV) and L_ 2 (As = 0.07 meV); (b) SL (A s = 
4.48 meV) and SL-i (As = 4.09 meV). 

[T2l The wave function of the ground state becomes 

|*SO£_ 1 ) -> II ( aC lx + bc lx + C 4,X + dc l,i) 47 ) 
X 

XC 7,X c lx\°) 

with 

|a| 2 + |6| 2 + |c| 2 + M| 2 -l. (148) 

We show in Fig. [12] how the different polarizations 
L z ,P z ,O z and O x change with bias in SX_i when E11 
is increased. The other polarizations are zero. The 
spin(layer) polarization S z (L z ) increases (decreases) with 
with E|| until > 0.2 mV/nm where it remains con- 
stant. The orbital coherence has not yet saturated at 
En. The ground state above this critical electric field is 
represented in the inset at the top-right corner of Fig. [12] 
It is interesting that the spin polarization can be varied 
in this phase by an external electric field. 

The state |^sol_i) has one gapless Goldstone mode 
and its dispersion is anisotropic in wave- vector space (as 
is the dispersion of the other modes). By contrast, if we 
apply E|| to a phase O v , the U(l) symmetry of the dipoles 
in the x — y plane is broken. The orbital pseudospin O 
is then forced to align with E|| and the Goldstone mode 
is gapped 2 ^ . 
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The absorption in l^soi-i) shows two peaks as in 
I^SL-j) • The first, low -energy, peak is shown in Fig. 
Il3l for two different orientations of the electromagnetic 
wave polarization. Clearly, the absorption is anisotropic 
in phase SOL—i by contrast to all the other uniform 
states that we studied before. 
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FIG. 12: (Color online) Variation of the spin and layer polar- 
ization and of the orbital coherence with an applied in-plane 
electric field in phase SOL-i. The bias As has been taken 
near the middle of the SL_i phase where (S z ) = (P z ) w 1. 
Parameters are B = 10 T and k = 5. 
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FIG. 13: (Color online) Electromagnetic absorption in the 
SOL-i phase for two different polarizations of the electro- 
magnetic field. 



VI. CONCLUSION 

In this work, we have derived the phase diagram of 
the C2DEG in a Bernal-stacked graphene bilayer. For 
the non-interacting Hamiltonian, we used a tight-binding 
model with the hopping terms 70 , 71, 74 and S and intro- 
duced a potential bias between the two layers given by 



A b ■ The Coulomb interaction was treated in the Hartree- 
Fock approximation. To reduce the complexity of the 
problem, we used an effective two-band model which de- 
scribed the low-energy behavior of the C2DEG and is 
valid for As << 71 ■ We also restricted the Hilbert space 
to the N — Landau level only and worked at zero tem- 
perature. Our method allows us to include both coherent 
and incoherent phases in the phase diagram. Indeed, we 
found phases with layer coherence at small bias, spin and 
layer coherence at intermediate bias and orbital coher- 
ence at large bias. The application of a parallel electric 
field, as we showed, can also lead to a state with orbital, 
layer and spin coherence. 

We have included in our analysis the hopping parame- 
ter 74 which is often neglected in theoretical calculations. 
In our calculations, we find that the phase diagram in 
sensitive to the precise value of this parameter. If the 
value of this parameter is modified, the phases that we 
have discussed are still present in the phase diagram but 
they occur at different bias. Moreover, other phases may 
appear. For example, with 74 = 0, the orbital-coherent 
phase becomes possible at v = — 1 and a new phase with 
valley and orbital coherence appears as discussed in Ref. 

HE 

We have written down the ground-state wave function 
for each phase in the global phase diagram of the C2DEG. 
We have also calculated for each phase the transport gap, 
the spin polarization, the collective mode dispersions and 
the electromagnetic absorption spectrum. The change in 
these properties from one phase to another should facili- 
tate their experimental identification. Strictly speaking, 
however, our results are only valid within the limits of 
validity of the approximations listed above. In partic- 
ular, we have neglected screening corrections which are 
known to reduce significantly the transport gaps. These 
corrections were considered for the incoherent phases in 
Ref. [l7l . In principle, these screening corrections should 
be smaller at larger magnetic field. The stability of the 
different phases that we found should also be studied by 
considering quantum and thermal fluctuations as well as 
disorder effects. 

More subtle corrections specific to graphene have also 
been considered by Shizuya 20 ~ ? - . The quantum fluctu- 
ations of the Dirac sea (the filled Landau levels from 
the valence band) have been shown to be sizable and to 
lead to corrections of the energy of the octet of states in 
TV = 1 — . According to Shizuya, the orbital degeneracy of 
the zero-energy levels is lifted by Coulombic vacuum fluc- 
tuations, leading to an appreciable shift and splitting of 
the n = and n = 1 levels and to a negative capacitance 
effect that blocks the rotation of the valley pseudospins. 
The negative capacitance effect appears when the full 
four bands of the tight-binding model are considered. 

A more complete calculation would include all these 
effects and allow a more direct comparison with the ex- 
perimental results. 
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